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Abstract 



We address the question whether we have to perform the analytic continuation to 
4 space-time dimensions before or after the perturbative expansion has been made. 
Using the simple model of large- A scalar field theory, we show how this affects quan- 
tities like the thermal mass and the pressure, why the two procedures give different 
results, and how this relates to the infrared behaviour of the theory. We conclude that 
the correct procedure is to take the limit d — > 4 before performing the perturbative 
expansion. Infrared divergences arise if one makes the perturbative expansion before 
going to the physical space-time dimension and they correspond to logarithms of the 
coupling in the exact (nonperturbative or resummed) expression. Finally, we outline 
how this may be relevant to the infrared problems encountered in thermal QCD. 



1 Introduction 



It has been known for a long time that beyond-leading-order calculations in mass- 
less thermal field theories exhibit infrared divergences fll], El- I* 1 the case °f abelian 
theories, systematic resummation techniques have been developed to cure infrared 
divergences. This leads to perturbative expansions which are nonanalytic in the cou- 
pling. In contrast to the case of abelian theories, no resummation method has so 
far managed to cure infrared divergences in non- abelian gauge theories like QCD 
PJ. In QED, infrared divergences mainly arise in Feynman diagrams with multiple 
insertions of the static photon self-energy. A resummation of all diagrams with an 
arbitrary number of such insertions (ring diagrams) can be shown to eliminate all 
infrared divergences. Consequently, quantities like the free energy can, in principle, 
be expanded to arbitrary order in powers of the coupling. In QCD however, since 
gluons couple to themselves, the behaviour of the theory is much more singular: di- 
agrams with static spatial gluon loops are increasingly infrared divergent when their 
order gets higher. This may signal the presence of a magnetic mass which would 
act as an infrared cut-off to regulate infrared divergent integrals. Although such a 
magnetic mass is potentially of order g 2 T since spacelike gluons remain massless to 
lowest order, it cannot be computed perturbatively. This is because beyond one-loop 
order, all the diagrams in a self-consistent expansion of the magnetic mass will be of 
the same order in the coupling 0. As a result, unresummed perturbation theory does 
not seem to work any more and all those diagrams that contribute to the same order 
in the coupling should somehow be resummed in a similar fashion to simpler theories 
like A0 4 or QED. Due to the complexity of QCD, this is unfortunately impossible in 
practice. 

If one uses dimensional regularisation, with d the dimension of space-time, the break- 
down of perturbation theory from some order signals itself through infrared l/(d — 4) 
poles in the perturbative series from a given order. It reflects the fact that the exact 
quantity one tries to compute is nonanalytic in the coupling when d — 4. In this work 
we therefore examine dimensional regularisation at finite temperature in more detail. 
Dimensional regularisation was devised as a method to separate the finite and diver- 
gent parts of Feynman diagrams arising beyond leading order in perturbation theory 
||. At finite temperature, there are infrared l/(d — 4) divergences in addition to the 
zero temperature ultraviolet ones so, if one use dimensional regularisation in the usual 
way, one has to make sure that these IR divergences somehow cancel each other in the 
expression for physical quantities. Consequently, one could wonder whether making 
an expansion in powers of g and d — > 4 commute. This is relevant because if we keep 
d — 4 nonzero, the perturbative expansion in powers of g actually involves gj{d — 4) 
and therefore is illegal since we are due to take d to 4 at the end of the calculation. 
As a result, an expansion in powers of g is is only valid if l/(d — 4) poles somehow 
cancel each other in the expression for physical quantities. A case where the previous 
limits do not commute is the following: consider 
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If we first expand A in power series of the coupling, its terms are divergent in the 
limit d — > 4 : 

C Cg 2 

However, in 4 dimensions, A = C/g 2 is perfectly finite and is nonanalytic in the 
coupling. 

This work investigates the consequences of this kind of problem in the case of large- 
N scalar theory, which will give us some insight about the way things might work 
in QCD. In next two sections, we derive and renormalise an exact formula for the 
thermal pressure in the case of large- N scalar field theory. In section 4, we study 
the singular structure of the expansion of the thermal mass and show that infrared 
divergences arise from a given order onwards in the renormalised perturbation series. 
We argue that those divergences crop up because the perturbative expansion has 
been performed before taking d to the physical dimension, in the same fashion as for 
the example given in We will see that the way to rescue the situation is to 

resumm the perturbative series and that we then obtain the correct expansion, ie the 
one we get if we go to 4 space-time dimensions first. Finally, we show this problem 
originates in the fact that the finite temperature contribution to the self-energy is 
actually nonanalytic in the coupling when d = 4 and involves terms proportional to 
g log g and half- integer powers of the coupling. 



2 Unrenormalised formula 



In this work, we use large- iV scalar field theory as a toy model to investigate the 
commutativity of the analytic continuation to four space-time dimensions and the 
perturbative expansion. Consider the theory described by the unrenormalised La- 
grangian density 

C(X,x) = hd^f - \m 2 ct> 2 - ^^0 4 (2.1) 
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where <fi is an iV-component field. We are interested in the limit iV 
fixed. The grand canonical partition function reads 
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where j3 = 1/T and r is some arbitrary time. The thermodynamic definition of the 
pressure is, as usual, 

p = TV' 1 log Z/3 (2.3) 

If we perform the change of variables <\> — > \ x l 2 <\) in the path integral, we can make 
the coupling dependence of the partition function more apparent M and obtain 



;2.4) 



z P {\) 


JV(f) exp 


i rr- 
X Jt 


- i(3 dx° J d d ~ 


x x £(1, x) 




/ V(p exp 


i rr- 
X Jt 


' ip dx° J d d ~ 


- l x £{0,x) 



If we define P = P(T, A) — P(T, 0) and use the translation invariance property that 
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it is straightforward to obtain from 



the unrenormahsed equation 
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Because it will simplify the equations in what follows we will hereafter use the notation 
AA =< A > — < A > free for any operator A. 



2.1 Calculation of AC 

The main analytic non-perturbative tools available here to compute AC are the 
Schwinger-Dyson equations which can be regarded as the quantum equations of mo- 
tion for Green's functions ||. These equations are the same as at zero temperature; 
the only difference lies in the boundary conditions : at finite temperature, the time 
periodicity (for bosons) or anti-periodicity (for fermions) of classical fields (ie the in- 
tegration variable in the path integral) provides Green's functions with temperature 
dependent boundary conditions. In large- TV scalar field theory, the self-energy is ex- 
actly calculable and it is therefore desirable to derive an expression for AC that only 
involves the propagator since, using the Dyson resummation, the latter can in turn 
be written in terms of the self-energy. To achieve this, we have to express the thermal 
four-point function in AC in terms of the propagator. This can be done using the 
Schwinger-Dyson equation for the propagator. To derive this equation, consider the 
generating functional of the thermal Green's functions : 
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We use the Keldysh variant of the real-time thermal field theory jl], |![], so the contour 
C runs along the real axis from — oo to oo, back to — oo and then down to — oo — i/3. 
According to the usual Feynman prescription, we have to add ie<j) 2 to the Lagrangian 
density in the generating functional fl2.1|) in order to ensure the convergence of the 
path integral for large <fi (in what follows, this prescription will be implicit). Then we 
have the following identity 
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which reads 
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By applying S/S(f>(y) to both sides of this equation and setting J to zero, we obtain 



(T c 0(y)(9 2 + m 2 )0(x)) + - 
and similarly 
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where T c and 5 C denote the time ordering and delta function along the integration 
contour in complex time. When x = y, the last two equations lead to 
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With the Keldysh contour, we get 



AC = - 
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where the subscript 12 refers to the element of the 2x2 thermal matrix propagator. 
This propagator has the structure 



D(g) = M(q°) D(g) M(q°) 



(2.9) 



with 
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where n(q°) is the Bose distribution (e 13 ^ — l) . This implies 
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In the case of large- scalar theory, the self-energy U(q, m, T) is real and independent 
of the momentum. In consequence, we have 
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Hence we find 
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where M, Nt and Mt are defined as : 

M (m 2 ) = - f % - = - / -^-27r5 + (q 2 -m 2 ) (2.14) 

V ; 2 J (2n) d q 2 - m 2 + ie 2 J (2vr) d W ) \ ) 

N T (m 2 ) = J 2vr S + (q 2 - m 2 ) n(q°) (2.15) 

M T {m 2 ) = M{m 2 ) + N T {m 2 ) (2.16) 

To make use of ( |2.13j ), we need to express all the bare parameters of the Lagrangian 
in terms of their renormalised counterparts. 



3 Renormalisation 

We choose a renormalisation scheme that renders AC as simple as possible. Since 
we are dealing with exact equations, it is suitable to define the renormalised mass 
through the Dyson equation, which is valid to all orders. At T = and in the on-shell 
renormalisation scheme, define the renormalised mass as 

m 2 R = m 2 + U(q 2 =m 2 R ,m,T = 0) (3.1) 

For large- iV scalar theory, the Dyson equation reads 

n(m, T = 0) = \M(m 2 R ) (3.2) 

so we have 

m \ = m 2 + \M{m 2 R ) (3.3) 
since Nt vanishes at T = 0. The integration in (|2.14|) gives 

"W> = T -^ffi < 2 <«> 

At nonzero temperature, it is convenient to define a thermal mass 

5m 2 R = X(M T (m 2 R + 5m 2 R ) - M(m 2 R )) (3.5) 

so ( |3.3| ) implies 

m R + 5m R = m 2 + XM T (m 2 R + 5m R ) (3.6) 

Comparing ( |3.6|) with (|3.3|) , we see that the definition (|3.5|) is a natural one. In what 
follows, we shall be interested in the case of massless theories (m = and therefore 
rriR = 0) because there one encounters the infrared problems of thermal perturbation 
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8m R = XM T (5m R ) (3.7) 
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Beside mass renormalisation, we also have to renormalise the coupling constant: be- 
cause A has dimension M 4 ~ d , we introduce a renormalisation scale // and define a 
dimensionless renormalised coupling by 



X = ^- d X R Z(X 
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where Z(Xr) is a combination of wave-function and vertex renormalisation factors. It 
is a function of Xr only because it is dimensionless. In the case of massive theories, 
one usually uses the on-shell renormalisation scheme 

Xr = A«Gu) = Z" 4 X + XX R M'(m 2 R ) (3.10) 

where the prime denotes differentiation of M with respect to its argument. However, 
if we consider massless theories, this definition is not useful since M'(0) is infrared 
divergent. Instead, define 



X R (fx) = /" 4 A + A A fl (/x) M'(/i 2 ) 
With (pD and (p|), this gives 
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In fact, the renormalisation scheme defined by ( |3.11|) is the modified minimal sub- 
traction scheme (MS) since ( [Up has the expansion 
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If we introduce the /3-function as usual as 
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and if we differentiate the left and right-hand sides of ( |3.9|) with respect to \i at fixed 
unrenormalised coupling, we get 
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From this equation and from ( |3.12| ), 
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We can also differentiate 
obtain 

-T- = d\ R 



with respect to \r at fixed /i and then use ( |3.15| ) to 
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The pressure being a physical quantity, it is the same before and after renormalisation: 
P R (T , \ R ) = P{T , A). We can therefore combine 0, d3A2|) and Q3J7D to turn (gj) 
into the renormalised formula 
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since 



This equation is of course consistent with its unrenormalised counterpart 
A = \r to lowest order and Pr = P. It also implies that the question whether 
we get the same answer for the thermal pressure if use perturbation theory before or 
after we take the d — > 4 limit reduces to how Sm R behaves depending on the chosen 
procedure. In the next section, we will find that if we make an expansion in power 
series of the coupling (ie use standard diagrammatic Feynman perturbation theory) 
in d space-time dimensions, not only we get uncancelled l/(d — 4) divergences from 
order X R onwards but even the finite second order term is incorrect when we compare 
it to the exact alternative expansion we obtain in 4 dimensions. 



4 The thermal mass 

We now study the structure of the mass shift when the coupling constant is chosen 
as expansion parameter either before or after we go to 4 dimensions. Our aim is 
to see what kind of expression we get in a perturbative calculation of the thermal 
mass and how resummation works to remove infrared l/(d — 4) terms and obtain the 
unrenormalised non perturbative equation (|3.7Q . 

4.1 Naive perturbative expansion 

In scalar field theories, it is known that the dominating infrared contributions to the 
self-energy come from the so-called ring diagrams and their extension called super- 
daisy or Hartree-Fock diagrams (see figure 1). 

First compute the value of a ring diagram with L external loops attached using the 
imaginary-time formalism : 

(L) AT 2L L ( ~ r d d - l p 1 r d d - l p 1 \ 

II [d,T) - — m T [-1) \2^J (2^-1^2 + ^+1+7 (2vr)^V L+ V 
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L loops attached 





Figure 1: Ring and Hartree-Fock Diagrams 

where u n denote the Matsubara frequencies and tut = Mr(O) 1 / 2 the one-loop thermal 
mass. According to the rules of dimensional regularisation, the last term in this 
expression, the static mode contribution, vanishes since the integrand does not contain 
any dimensionful quantity. However, the sum over L of the static mode contributions 
is 

V ^ L (-1) L ( J^J_ = ^ f 1 = U 2) 

2 T{ ' J (2 7 r) d - 1 p 2L + 2 " 2 J (27TY- 1 p 2 + m 2 T " 8tt 1 ' 

The perturbative expansion of the right hand side of (|4.2|) is inappropriate since it is 
nonanalytic in the coupling (it is proportional to A 3 / 2 ) and naive use of dimensional 
regularisation fails here. We can therefore conclude that static mode contributions 
only make sense when they have been resummed. 

We now turn to the non-static modes. In the spirit of dimensional regularisation, 
integrals are initially computed in a dimension where they converge and then analyt- 
ically continued to four space-time dimensions. In d dimensions, the non-static part 
of the self energy is 

nSK T) = XT m?(-l) L {2 *%T 3 T{L + 1 ~ ¥) C(2L + 3 - d) (4.3) 

(47f)~ l - 

where ( is the standard Riemann zeta-f unction. As usual, we get to lowest order 

XT 2 

*4 = -n£i.(4,T) = ^L (4.4) 

Since ((z) is analytic in the whole complex plane except for z = 1 where it has a 
simple pole, njf] (d,T) is finite for all values of L when d — ► 4, except for L — 1. In 
the latter case, we have 

n « ( d , T) = -^i(2-Ilog47r-logT)+ 0(d - 4) (4.5) 

n - sA ' ; 167r 2 (d-4) 16vr 2 V2 2 & B J v ; K J 

It means we still get a X 2 /(d — 4) pole even if we have resummed all the static 
contributions. 



Since it is known that resumming all ring diagrams provides a finite result ||, the 
resummation of all the Il( L ) (d, T) for L > 1 should cancel out the — 4) pole 
coming from (d,T) : 



AT 



(4tt) - 



/ m T \ 2L (-l) L r(L 



L=2 j n=X 

If we use the Taylor expansion 
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we find, by summing over n, 
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By expanding this expression around d — 4, we obtain 
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As expected, resummation indeed removes the pole: it is cancelled out by the rest of 
the perturbative series. Hence, the resummed result is merely the one-loop self-energy 
diagram in which the free propagator has been replaced by the first-order dressed one: 



J2 n (i) {d, T) = A M T (m 2 T ) + 0(d - 4) 



(4.11) 



L=0 



Because it comes from the Matsubara sum, it is not clear whether the A 2 / (d — 4) 
divergence was an infrared or an ultraviolet one. In order to find out, it is interesting 
to study nW(d, T) in the real-time formalism: this is done in appendix A.l. 
In this section, we have so far only been dealing with quantities involving the unrenor- 
malised coupling so we still need to renormalise (|3.7|) to make sure we have properly 
taken care of UV divergences before jumping to conclusions about possible remaining 
divergences. Expressing bare quantities in terms of renormalised ones, ( p.7|) becomes 
the following integral equation : 



5m 2 R = /i 4 d X R 



M(fj, 2 ,Sm 2 R ) + N T (Sm 2 R ) 



(4.12) 



n 



where M(fi 2 ,8m 2 R ) = M(5m 2 R ) — 8m 2 R M'(fi 2 ) is finite when d — 4. This equation 
is exact and therefore corresponds to the full summation of all the Hartree-Fock dia- 
grams. By studying the ring diagram contributions to the thermal mass in the ITF 
and RTF, we have seen that an expansion of NT(8m 2 R ) in power series of 8m 2 R gives 
rise to infrared divergences of the form Sm R / (d — 4) although itself is perfectly 
finite when d — 4. This is because Nt is nonanalytic in 8m R at 8m\ = (namely, it 
involves terms which are proportional to 8m R log 8m r) and so the Feynman perturba- 
tive expansion in powers of X R is illegal since 8m R ~ X R to lowest order. In the next 
section, we investigate the exact form of Nt and its nonanalyticity in the coupling. 
This will lead us to the conclusion that the 8m 2 R j [d — 4) term in the perturbative 
expansion of Nt is effectively replaced with 8m\ log(8m 2 R ) in the corresponding exact 
expression (see ( 4.18 )). In section 4.3, we shall see that the divergent term in a per- 
turbative expansion of NT(8m 2 R ) gives rise to infrared divergences in the renormalised 
perturbative series from order X R onwards. 



4.2 Alternative expansion of TV 
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In order to study its dependence in the the coupling, we will now expand Nt in 
another and more sensible way than the naive perturbative expansion ||: 
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where y = (3m. We evaluate the sum using the Mellin summation formula (see 
appendix A. 2) 

00 1 r C+ioo 

Ef( k ) = W- / dz((z)cf>(z) (4.14) 

where <j>{z) is the Mellin transform of f(k) and C is such that all the poles of the 
integrand are to the left of the vertical integration line. We obtain 
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2i7r(47r) — T( 



C—ioo 



with C > d — 2. We can evaluate the integral over u and the remaining integral 
becomes 
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The integrand has single poles at z = 0, 1 and 2; = d — 2k where = 1,2,3,.... The 
contour of the z integration in the integral (|4.15|) may be closed in the left half plane 
since the integrand vanishes asymptotically on the arc at infinity. We now have to 
compute all the residues of the poles of the integrand. The contribution to iVV from 
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the pole at z = 2 — d corresponds to the first order thermal mass m\. From the pole 
at z — 1, we get — Tm/(87r) which corresponds to the sum of all the static modes we 
had in the ITF (see ( |4.2[) ) and is nonanalytic in the coupling (mj ~ A]/ 2 )- The sum 
of the contributions to N? from the poles at z = and z = d — 4 is 
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and this is nonanalytic in the coupling as well due to the term in logm. Finally, all 
the remaining poles have residues of the same form and the complete expression for 
Nt reads (when d = 4) 
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From this expression, we see that NT^m 2 ^) contains a term involving logA# and 
another one proportional to \]( 2 . The fact that an expansion of Nr r {8m 2 R ) in power 
series of the coupling is illegal is quite clear from ( 4.13|) : if we try to expand the 
integrand in power series of m 2 , the expansion parameter is actually m 2 /p 2 and this 
obviously does not make sense since the momentum integration runs down to 0. If we 
do it anyway, we obtain what we got in perturbation theory and we get an isolated 
l/(d — 4) term at order X 2 ^. Moreover, we also have 
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So we find 
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One can solve this equation iteratively to any desired accuracy in Xr and the first few 
terms in the expansion are [0] 



Sm 2 R 

J^2 



Xr 
24 



A 



3/2 



16ttV6 
l - 27 - 2 lo, 



+ 3 - 7 - log 



\ 2 



A-nTJ 384tt 2 



A 



5/2 

R 



4ttT/ i0247r 3 J2/3 



+ 0(X%) 



(4.21) 



It is evident that Sm 2 R does not have an expansion in powers of Xr, but rather in 
powers of X l R 2 . 



4.3 Renormalised perturbative series 



We express 5m R in terms of the renormalised coupling and the renormalisation scale. 
From Q and (jUll) 
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so diagrams corresponding to different fixed orders in A will yield contributions both 
to the same order in X R and to higher orders. Evidently, to lowest order, we have 
A = X R so Sm R = X R T 2 /24. To order A^, both the one and two- loop diagrams 
contribute and the l/(d — 4) divergence coming from II^ 1 ) exactly cancels the one that 
results from the renormalisation of the coupling in the one-loop diagram n(°). We are 
therefore left with a finite expression up to order A^: 
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Having obtained an exact expression for Sm 2 R in terms of Xr, we investigate the 
problems associated with first expanding it as a perturbation series in Xr. Comparing 
with ( |4.21| ), we see that there are missing terms. This is because we also need to 
include the contribution from the non-vanishing sum of all the static modes : this 

1/2 

contribution, which also yields all the terms proportional to an odd power of X R in 
Q4.21D , is missed out in the perturbative expansion and only appears when we resumm 
the series. 

After tedious calculations, divergences cancel out again at order A^ if we take into 
account all the diagrams that contribute to the thermal mass at this order (they 
are represented in figure 2): diagram D in figure 2, which is not a ring diagram, 
rescues the situation by cancelling divergences coming from diagrams A and B and 
the renormalisation of the bare coupling in these diagrams. It means that, if we only 
consider ring diagrams, we get uncancelled divergences from order X 3 R onwards in the 
renormalised perturbation series. From order A^. onwards, it is not possible to get 
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Figure 2: Diagrams contributing to Sm R at order A^ 
rid of divergences any more even if we sum up all diagrams that contribute to the 



thermal mass at this order : worse and worse divergences crop up in the renormalised 
perturbative series when going to higher orders. For instance, at order Afj, we have 



Sm 2 R 



o(x%) 49152vr 6 (d-4) 



A;// ~ • 0{{d- 4) -) (4.24) 



In conclusion, we have seen in this section that the nonanalyticity of the thermal 
mass in the coupling causes the failure of the corresponding Feynman perturbative 
expansion. This failure has two consequences : firstly, infrared divergences appear in 
the renormalised series from order Afj and, secondly, the term proportional to A 1 / 2 
and its contributions to all order in A# are missed out. 



5 Conclusion 



The equation ( |3.18|) derived for the pressure is exact, or non perturbative. Its right- 
hand side involves Sm^ and we have seen in the previous section that we get uncan- 
celled l/(d — 4) infrared singularities if we expand it in powers of Xr before we take 
the d — > 4 limit. On the other hand, if we take the limit d — > 4 on the exact expression 
( [4.13[ ), then we get an expression which cannot be expanded in power series of the 
coupling any more since a coefficient in the expansion would be infinite. It tells us 
that we have to use the more involved expansion outlined in section 4. In fact, before 
we take the limit d — > 4, nothing prevents us from expanding in power series of the 
coupling even if this expansion is illegal in four space-time dimensions. We could 
therefore conclude by saying that the correct procedure is to work in four space-time 
dimensions before expanding the physical quantity under consideration in a suitable 
way. In this paper, we studied both the perturbative and the exact expression of the 
finite temperature contribution to the thermal mass and we saw that the infrared di- 
vergent A 2 / (d — 4) term in the unresummed perturbative series is effectively replaced 
with A log A in the corresponding exact expansion (see section 4.2). The same kind 
of mechanism may be related to the severe infrared problems encountered in thermal 
QCD. In the case of the QCD thermal pressure, it is known that diagrams with mag- 
netostatic gluon loops become divergent from order g e in perturbation theory. From 
what we saw in this work in the case of large- N scalar theory, it seems sensible to 
assume that this g 6 / (d — 4) divergence in the perturbative QCD pressure corresponds 
to a term involving g 6 log g in the exact expression for the pressure. This will be the 
subject of subsequent work. 
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A Appendix 



A.l Ring diagrams in the real-time formalism 

In the real-time formalism, the self-energy also has a matrix structure and the matrix 
element relevant to the evaluation of (d, T) is 



n$(d,T) = -i±M> 
Using ( j2.9Q , we get 

dIM-dUp)dIM 
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The finite temperature part of the latter expression can be written as 

d f I 1 
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m 2 =0 



which is the derivative of the one-loop diagram for the thermal mass. We therefore 
obtain 



d d p 



D\M ~ D 2 12 {p)DlM =-i7r(d-3) 



Finally, if we use the representation 
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and the relation 



C(l-*) = 2(27r)-*cos( T ) T(z)C(z) 
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(A.4) 



(A.5) 



(A.6) 



we obtain the same result ( [4.5[ ) as with the imaginary-time formalism. Since the mo- 
mentum integral in the right-hand side of (|A.4j ) is infrared divergent in four space-time 
dimensions, it tells us that the corresponding l/(d — 4) divergence really is an infrared 
one. Similarly, we could also rederive the general expression for a ring diagram with 
an arbitrary number of external loops and check that it agrees with what we found 
using the imaginary-time formalism : one just has to take higher derivatives with 
respect to m 2 in ( |A.3| ). 



A.2 The Mellin summation formula 

The Mellin transform of a function f(x) is usually defined by 

<j>( z ) = / x z ~ 1 f(x)dx (A.7) 
Jo 

1 I 



The transform 4>(z) exists if the integral 



/ \f{x)\x k - l dx 
Jo 



(A.8) 



is bounded for some k > 0, in which case the inverse transform is 



C+ioo 



C'—ioo 



x z <fi(z) dz 



(A.9) 



where C > k. Together with the standard definition of the Riemann Zeta-function, 
this gives the summation formula 
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